State space structure and entanglement of rotationally invariant spin systems 
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We investigate the structure of SO(3)-invariant quantum systems which are composed of two 
particles with spins ji and J2 . The states of the composite spin system are represented by means 
of two complete sets of rotationally invariant operators, namely by the projections Pj onto the 
eigenspaces of the total angular momentum J, and by certain invariant operators Qk which are built 
out of spherical tensor operators of rank K. It is shown that these representations are connected 
by an orthogonal matrix whose elements are expressible in terms of Wigner's 6-j symbols. The 
operation of the partial time reversal of the combined spin system is demonstrated to be diagonal 
in the Qx-representation. These results are employed to obtain a complete characterization of spin 
systems with j\ — 1 and arbitrary j'2 > 1. We prove that the Peres-Horodecki criterion of positive 
partial transposition (PPT) is necessary and sufficient for separability if j% is an integer, while for 
half-integer spins j% there always exist entangled PPT states (bound entanglement). We construct 
an optimal entanglement witness for the case of half-integer spins and design a protocol for the 
detection of entangled PPT states through measurements of the total angular momentum. 

PACS numbers: 03.67.Mn,03.65.Ud,03.65.Yz 



I. INTRODUCTION 

Entanglement is a basic feature of composite quantum 
systems connected to the tensor product structure of the 
underlying Hilbert space of states. A mixed state of a 
bipartite quantum system described by some density ma- 
trix p is said to be entangled or inseparable if p cannot be 
written as a convex linear combination of product states. 
Otherwise it is called classically correlated or separable 
Q. The properties of entangled states are responsible 
for many of the fascinating and curious aspects of the 
quantum world and lie at the core of many proposed ap- 
plications in quantum information processing 

The general characterization and quantification of en- 
tanglement in mixed quantum states is a highly non- 
trivial problem. It is even very difficult in general to for- 
mulate simple operational criteria which allow a unique 
identification of all separable states of a given composite 
system. There do exist, however, many necessary sepa- 
rability criteria 0, 0, S, 0, E3, E3, E3, E3| - A simple 
and, in fact, very strong criterion is the Peres-Horodecki 
criterion 0, 0] which states that a necessary condition 
for a given density matrix p to be separable is that it 
has a positive partial transposition (PPT states). It is 
known that this criterion is necessary and sufficient for 
certain low-dimensional systems, while it is only neces- 
sary in higher dimensions Q. 

The analysis of the entanglement structure is greatly 
facilitated through the introduction of symmetries, i. e., 
if one restricts to those states of the composite system 
which are invariant under certain groups of symmetry 
transformations. Important examples in this context are 
the manifolds of the Werner states 0], of the isotropic 
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states 0, Q] and of the orthogonal states ^^|. Here, 
we investigate entanglement under the symmetry group 
SO (3) of proper three-dimensional rotations of the coor- 
dinate axes. More precisely, we consider the problem of 
mixed state entanglement of systems which are composed 
of two particles with spins j\ and j2 , and which are invari- 
ant under product representations of the group SO(3) or, 
equivalently, of the covering group SU(2). A basic tool 
of our analysis is the work of Vollbrecht and Werner 
which provides a general scheme for the treatment of en- 
tanglement under given symmetry groups. 

Mixed SO(3)-invariant states of composite systems 
arise, for example, from the interaction of open systems 
[lf3 | with isotropic environments Their analysis is 

of great importance and leads to many applications. As 
examples we mention investigations on the connection 
between quantum phase transitions and the behaviour 
of entanglement measures (see, e. g., [011111]), the analy- 
sis of entanglement of SU(2)-invariant multiphoton states 
generated by parametric down-conversion [2fj|, and stud- 
ies of the entanglement of formation [2J. The technique 
of this paper could also be relevant for the characteri- 
zation of quantum correlations in Fermionic or Bosonic 
systems developed recently USES]- 

The Hilbert space of a system which is composed of 
two particles with spins j\ and ji is given by the tensor 
product C Nl ® C^ 2 , where N x = 2ji + 1 and N 2 = 2j 2 + 1 
are the dimensions of the local spin spaces. We call such 
a system an Ni <g> N 2 system. Throughout the paper we 
will assume that j\ < j 2 , i- e., N\ < N 2 . 

According to the Peres-Horodecki criterion 0, Q the 
cases of 2 ® 2 and 2 ® 3 systems are trivial: It is known 
that in these cases the PPT criterion is necessary and 
sufficient for all states, i. e., even for states which are not 
invariant under rotations. Schliemann |2J] has shown 
recently that the PPT criterion is also necessary and suf- 
ficient for SO(3)-invariant 2®N 2 systems with arbitrary 
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N2 . The case of 3 ® 3 systems has been treated by Voll- 
brecht and Werner [Ta j, who proved that the PPT cri- 
terion is again necessary and sufficient for separability. 
For 4 ® 4 systems a qualitatively new situation arises: It 
has been demonstrated in (2^1 that the PPT criterion is 
not sufficient and that the entangled PPT states form 
a three-dimensional manifold which is isomorphic to a 
prism. In the present work we investigate the important 
special case of 3 (g> N2 systems with arbitrary iVjj . 

The method developed in [5I| enables the treatment 
of the case of equal spins ji = J2. In this paper we 
extend this method to arbitrary spins ji and 32. For 
the analysis of entanglement under SO(3)-symmetry it 
is advantageous to replace the transposition used in the 
PPT criterion by another unitarily equivalent operation, 
namely by the antiunitary transformation of the time re- 
versal. The reason for this fact is that the operation of 
the time reversal of states commutes with the represen- 
tations of the rotation group. 

There are two natural representations of rotationally 
invariant states. The first one uses the fact that any in- 
variant state can be written as a unique convex linear 
combination of the projections Pj onto the eigenspaces 
of the total angular momentum J of the composite spin 
system. The advantage of this representation is that it 
leads to very simple conditions expressing the positivity 
and the normalization of physical states. However, the 
set of the PPT states is most easily determined in another 
representation which employs the irreducible spherical 
tensor operators of spin-j particles. We will construct 
a complete system of invariant operators Qk which are 
built out of the spherical tensors of rank K. Any in- 
variant state of the composite spin system can then be 
written as a unique linear combination of the Qk- The 
introduction of the invariant operators Qk generalizes 
the ideas of Schliemann (2J, |23], who has developed a 
representation of SU(2)-invariant states by means of spin- 
spin correlators and has formulated various separability 
conditions and sum rules in terms of these correlators. 

The paper is organized as follows. The representations 
of SO(3)-invariant states in terms of the invariant op- 
erators Pj and Qk are constructed in Sec. |n] We also 
derive in this section the linear transformation which con- 
nects these representations and show that it is given by 
an orthogonal matrix whose elements are determined by 
Wigner's 6-j symbols. The behaviour of states under 
partial time reversal and the construction of the set of 
the invariant separable states are discussed in Sec. lIIII 

The general theory is then applied in Sec. IIVI to the 
case of 3 ® N2 systems with arbitrary N2 . We prove that 
the PPT criterion represents a necessary and sufficient 
separability condition for 3 ® N2 systems if and only if 
N2 is odd. Thus, for integer spins 32 all PPT states 
are separable, while for half-integer spins ji there always 
exist entangled PPT states. This fact has already been 
conjectured by Hendriks {13] on the basis of a detailed nu- 
merical investigation. We also show that for half-integer 
32 the boundary of the separability region is curved. Fi- 



nally, SeclVlcontains a discussion of the results and some 
conclusions. In particular, we construct an optimal en- 
tanglement witness for the case of half-integer spins and 
exploit this witness to design a protocol which allows the 
detection of entangled PPT states through measurements 
of the total angular momentum. 



II. REPRESENTATIONS OF SO(3)-INVARIANT 
STATES 

We consider two particles with spins j\ and 32 and cor- 
responding angular momentum operators j^- 1 ' and . 
The Hilbert space C Nl of the first particle is spanned by 
the common eigenstates |ji,mi) of the square of j'W and 

of jz \ where N\ = 2ji + 1 and mi = — ji, . . . , +ji- Cor- 
respondingly, the Hilbert space C^ 2 of the second particle 
is spanned by the eigenstates Ijz^mz), where N2 = 2j2 + l 
and m 2 = -32, • • ■ , +ja- 

The Hilbert space of the total system composed of both 
particles is given by the tensor product C Nl ® C* 2 . The 
angular momentum operator of the composite system is 
defined by: 



J = j-(l) 



I + 12) j 



(2) 



(2.1) 



where I denotes the unit matrix. A state of the composite 
system is described by a density matrix on the product 
space, i. e., by a positive operator p on C^ 1 ® C^ 2 with 
unit trace: p > 0, trp = 1. 

The irreducible unitary representation of the group 
SO (3) of proper rotations R on the state space of a parti- 
cle with spin 3 will be denoted by D^\R). The transfor- 
mation of the states of the composite spin system is then 
given by the product representation D^{R)<S>D^ 2 '{R), 
A state p of the combined system is said to be rotation- 
ally invariant or SO(3)-invariant if the relation 



D Ul >(R)(g> D {J2 >(R) p D^ 1 '(R)<g>D^ 2 >(R) 



-){.i-2 



)(j2 



P 



holds true for all proper rotations R € 5*0(3). 

We shall use two different representations of rotation- 
ally invariant states. The first one employs the projection 
operators 



+J 



Pj= Y, \JM){JM\, 



(2.2) 



M=-J 



where | JM) denotes the common eigenstate of the square 
of the total angular momentum J and of its z-component 
J z , i. e., we have J 2 1 JM) = J(J+1)\JM) and J Z \JM) = 
M\JM). The operator Pj projects onto the manifold 
which is spanned by the eigenstates belonging to a fixed 
value J of the total angular momentum. According to 
the triangular inequality J takes on iVi different values 
which may be integer or half-integer valued: 



J = 32 - 31, h - ji + 1, ■ ■ ■ , h + h- 



(2-3) 
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It follows from Schur's lemma that any invariant state p 
can be written as a linear combination of the Pj : 

P=^=Y^^Pj. (2.4) 

VT^Thy V2JT1 y ' 

Here, the aj are real parameters and we have intro- 
duced convenient normalization factors of y/WjNz and 
\]2J + 1. In order for Eq. 112. 4|l to represent a true den- 
sity matrix the aj must of course be positive and nor- 
malized appropriately: 

aj > 0, (2.5) 
x I2J + 1 . . 

^ = eva^= l (2 - 6) 

Any invariant state p is thus uniquely characterized by 
a real vector a. in an A^-dimensional parameter space 
IR^ 1 which will be referred to as a-space. The condi- 
tions of the positivity and of the normalization of p are 
expressed by the relations 112. 5 J) and H2.fi Jl . We denote 
the set of all vectors a whose components aj satisfy 
these relations by S a . Being isomorphic to the set of 
invariant states, S a is of course a convex set. We infer 
from Eqs. 112. 5 j) and 112. fij) that S a represents an (N± — 1)- 
dimensional simplex. 

A useful alternative representation of the invariant 
states is obtained by use of a complete system of irre- 
ducible spherical tensor operators (see, e. g. (2M l29l|l. 
The tensor operators which act on the state space of the 
particle with spin j% are written as , where i = 1,2. 
The index Ki = 0, 1, . . . , 2ji denotes the rank of the ten- 
sor operator. For a given rank Ki the index qi takes on 
the values qi = —Ki, —Ki + 1, . . . , +Ki. We thus have 
(2Ki + 1) tensor operators of rank Ki which trans- 

form under rotations according to an irreducible repre- 
sentation of the rotation group. The explicit definitions 
of the tensors and a brief summary of their properties are 
given in Appendix lAl 

Using the tensor operators one defines Hermitian op- 
erators Qk acting on the state space of the composite 
spin system: 

+K 

®K= E T k]® T k^ (2-7) 

q=-K 

where the index K takes on N\ different integer values: 

K = 0,l,...,2j 1 . (2.8) 

It follows from the transformation properties of the ten- 
sor operators that all Qk are invariant under rotations. 
For instance, the operator Qq is proportional to the iden- 
tity, Qo — -j=—I ® /, while Qi is proportional to the 

invariant scalar product j W . j( 2 ) of the spin vectors. 

The Qk defined by Eq. I|2.7I) form a complete system 
of operators. This means that any rotationally invariant 



Hermitian operator can be represented as a unique linear 
combination of the Qk in a way analogous to Eq. I|2.4I) : 

Here, we have again introduced appropriate normaliza- 
tion factors and real parameters /3k which form a vector 
f3 in an A^-dimensional parameter space IR^ 1 referred 
to as /3-space. The operators Qk satisfy ^{QkQk'} = 
(2K + 1)5 kk'- This fact follows directly from the or- 
thogonality relation ijAlll for the spherical tensors. The 
Qk for K ^ are therefore traceless which leads to the 
normalization condition 

trp = A, = 1. (2.10) 

The sets {Pj} and {Qk} represent complete systems 
of invariant operators. The corresponding parameter vec- 
tors a and (3 must therefore be related by a linear trans- 
formation R^ 1 1 ^ K Nl . We write 

(3 = La, (2.11) 

where L is an (JVi x N±) matrix. To find the elements of 
this matrix we use Eqs. 12.41 and 12.91 to get 

E -7%=PJ = E -J^=Qk. (2.12) 
Y V2J+1 Y V2^ + l 

Multiplying this equation by Qk 1 and taking the trace 
we find that the elements of L are given by 

L KJ = [(2K + 1)(2J + l)]- 1/2 tY{Q K Pj}. (2.13) 

This can be expressed as 

L K j = y/(2K+l)(2J+l)(-iy^' +J I I ]{ J K ]- 

(2.14) 

The curly brackets denote a 6-j symbol introduced by 
Wigner JsTj] into the quantum theory of angular mo- 
mentum. A proof of the relation 112.1411 is given in Ap- 
pendix [Bj The 6-j symbols are scalar quantities which 
are defined through invariant sums over products of 
Clebsch-Gordan coefficients. They describe the trans- 
formation between different coupling schemes for the ad- 
dition of three angular momenta j23|. Their properties 
have been studied in great detail and many closed for- 
mulae, recursion relations and sum rules are known. In 
particular, it follows from the sum rules that L represents 
an orthogonal (Ni x Ni) matrix. 

The above results lead to the conclusion that the set 
of SO(3)-invariant states is represented in /3-space by the 
set 

S 13 = LS a . (2.15) 

The set S 13 is again an [N\ — l)-dimensional simplex 
which may be constructed by determining the images of 
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the extreme points of S a under the orthogonal transfor- 
mation L. 

The introduction of two parameter spaces is motivated 
by the following observations. On the one hand, the set 
of states is most easily constructed as a subset in a- 
space. This is due to the fact that the representation 
of Eq. H2.4J1 corresponds to the spectral decomposition 
of p and, therefore, the requirement of the positivity of 
p immediately leads to the simple condition 112.511 . On 
the other hand, the representation IJ2.9I) of states in (3- 
space is much more suitable for the construction of the 
set of separable states, which is due to the fact that the 
operation of the partial time reversal is diagonal in the 
Q_r- -represent ation . 



III. INVARIANT SEPARABLE STATES 

A state p of the composite spin system is said to be 
separable if it is possible to write this state as a convex 
linear combination of product states: 



P : 



5> 



„(!) 



(2) 



A 4 >0, 53 Ai = 1 ' 



where the p| and p\ are normalized states of the first 
and of the second spin, respectively Q]. It is clear that 
the set in /3-space which represents the invariant and sep- 
arable states is a convex subset of S@ . This subset will 
be denoted by S@ cp . 

Following the work of Vollbrecht and Werner ^| we 
define a projection super-operator (SO (3) twirling) by 
means of 



ITp = dR U{R)pU{R)\ 



(3.2) 



where U(R) = D^(R)2>D^ 2 \R) and the integration is 
extended over all group elements R G 50(3). The twirl 
operation maps any state p of the composite spin system 
to an S0(3)-invariant state lip. Moreover, if p is separa- 
ble then also lip is separable. In terms of the invariant 
operators Pj or Qk the action of the twirl operation may 
be expressed by 



y tr{Pjp} = y HQkP} 

y J —f 2 J + 1 J ^ 2K + 1 ^ 



(3.3) 



It is known that any invariant separable state is a con- 
vex linear combination of ITprojections of pure product 
states. Given a pure product state 



p = I^ ( V 2) ><¥> ( V 2) I 



(3.4) 



Eq. i'.'i.'M shows that the corresponding parameters aj 
and Pk of its projection lip are given by 



a j 



2J+ 1 



(^V 2) I^MV 2) ), (3.5) 



fix = \/#fl^ ( V 2) |Q*MV 2) >. (3.6) 



We introduce into Eq. 13. (it the definition 112. 7J1 of the 

Qk and define the functions 



(2)1 



(3.7) 



= \/l^E^ (1) i^M 1) )^ (2) i4 2 i t M 2) >. 

^ q=-K 

Let us further define W 13 as the range of the parameter 
vector (3 whose components are given by these functions, 
where \ip (1) ) G C Nl and \ip (2) ) G C^ 2 run independently 
over all normalized states of the first and of the second 
spin, respectively: 

= {0\0K=0K[<P m ,<P (2) ], ||^ 2) || = l}. (3.8) 

The set of separable states is then equal to the convex 
hull of W p : 



Sf ep = hull (W?) . 



(3.9) 



This means that S^ ep is equal to the smallest convex set 
which contains W° . 

Within this formulation the problem of constructing 
S[~ ep reduces to the determination of the convex hull of 

the range of the functions (3 k- Even for the present case 
of a highly symmetric state space this is, in general, an 
extremely difficult task. A strong necessary condition for 
separability is the Peres-Horodecki criterion 0,0. Ac- 
cording to this criterion a necessary condition for a given 
state p to be separable is that its partial transposition is 
a positive operator: T 2 p = {I®T)p > 0. Here, TB = B T 
denotes the transposition of the operator B on C^ 2 which 
is defined in terms of the basis states of the second spin 
by means of (72, m 2 1 B T \j 2 ,m' 2 ) = (j 2 , m' 2 \B\j 2l m 2 ). The 
partial transposition T 2 is then defined by Ti{A % E) = 
A(g>B T . 

The operation of taking the partial transposition de- 
stroys the rotational invariance of states, i. e., if p is 
invariant under rotations the partially transposed state 
T 2 p is generally not SO(3)-invariant. However, there ex- 
ist another operation which is unitarily equivalent to T 2 
and which does map rotationally invariant operators to 
rotationally invariant operators. This operation will be 
denoted by z? 2 = It involves the antiunitary time re- 

versal transformation d of the second spin and will there- 
fore be referred to as partial time reversal. 

According to Wigner's representation theorem [j^l the 
action of the time reversal transformation d on an oper- 
ator B can be expressed as: 



■SB = VB T V^ = tB^< 



(3.10) 



In the first expression T denotes again the transposition 
and V is a unitary matrix which represents a rotation 
of the coordinate system about the y-axis by the angle 
7r. In the second expression of Eq. H3.1DJ1 t denotes the 
operator r = Vtq which is composed of the 7r-rotation V 
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and of the operator tq of the complex conjugation. The 
operator r is antiunitary and satisfies 



(3.11) 



■d is a positive but not completely positive map. It is 
unitarily equivalent to the transposition T and, hence, 
the Peres-Horodecki criterion can be expressed by 



$ 2 p = (I <8>tf)p > 0. 



(3.12) 



A great advantage of the representation of states in 
/3-space is that the operators Qk have a very simple 
behaviour under the map $2. Namely, as is shown in 
Appendix El they are eigenoperators of the partial time 
reversal: $ 2 Qk = {-1) K Qk- In (3 -space the map $ 2 
therefore induces a reflection of the coordinate axes cor- 
responding to the odd values of K: 



$2: Pk^ (-1) K P 



K ■ 



(3.13) 



We thus get the image $ 2 S 13 of S@ by reversing the signs 
of the odd coordinates. 

We define Sp pt as the set of states which are positive 
under d 2 or, equivalently, under T 2 (PPT states). This 
set is equal to the intersection of S@ with its image , d 2 S 13 . 
According to the Peres-Horodecki criterion the set of sep- 
arable states is a subset of the set of PPT states. Hence, 
we have 



(3.14) 



We note three properties which turn out to be useful 
in the construction of the set of separable states. 

(1) The functions defined by Eq. I|3.7ll are invariant 
under simultaneous rotations of the input arguments: 



Pk[d {31 \r)v [1 \d^\r)v {2) 



(3.15) 



This property is an immediate consequence of the rota- 
tional invariance of the operators Qk ■ 

(2) The range W 13 defined in Eq. 113. 8jl is obviously 
invariant under the partial time reversal i?2 • This means 
that /3 G W 13 implies tf 2 /3 € W 3 . 

(3) There exist two distinguished separable states. 
These are the state given by the parameter vector a with 
components 



N X N 2 



2c7rr 



1 



>J,J a 



J n 



h+j 2 , (3.16) 



and the partially time reversed state given by a' = r d 2 a. 
To proof this statement we consider a pure product state 
p of the form of Eq. I|3.4j) with {ip^) = \ji,+jx) and 



|^)) = 

\J = Jn 



|.72 , +J2) - We then have the obvious relation 



J 11 



\(p^(p^) and, hence, 



(3.17) 



Equation Ij3.5l) then immediately leads to Eq. H3.16J) . 
This means that the pure product state p is mapped 



under the twirl operation to the separable state lip = 
2J 1 +1 -P/ max corresponding to the maximal value of the 
total angular momentum Jmax- It follows from point (2) 
that also the partially time reversed state is separable. 

The point a given by Eq. H3.16J1 is an extreme point 
of the simplex S a and its image ex 1 is an extreme point 
of , d 2 S a . Thus, a. and a! are extreme points of Sp pt . 
It follows that the corresponding points (3 = La and 
/3' = La' in /3-space belong to W 13 and represent extreme 
points of S'ppf 

As an illustration of the above analysis consider a 
2 (g> N 2 system for which ji = \ and j o is arbitrary. As 
has been demonstrated by Schliemann j2J| the PPT cri- 
terion is a necessary and sufficient separability condition 
in this case. Within the present formulation this state- 
ment can be proven as follows. We first note that the 
index K takes on the two values K = 0, 1 such that (3 is 
a two-dimensional vector. Because of the normalization 
condition (12 . lOf) we only need a single parameter /3\ to 
characterize uniquely an invariant state of a 2 ® N 2 sys- 
tem. It follows that S@ can be represented by an interval 
of the /?i-axis, and Sp t by a sub-interval of this inter- 
val. Since an interval has exactly two extreme points (its 
endpoints) we conclude with the help of point (3) above 
that the extreme points of S^p,. belong to W 13 . By the 

relation l|3.9fl the sets 5p pt and S^ ep therefore coincide. 
This shows that the PPT criterion is indeed necessary 
and sufficient for separability. 



IV. 3 <g> N SYSTEMS 

Let us now consider the case ji = 1 (N% = 3) and j 2 
arbitrary, i. e. the case of 3(&N 2 systems. For convenience 
we write N = N 2 = 2j 2 + 1. Since J takes on the values 
J = j 2 — 1, j 2 and j 2 + 1, a is a three- vector 



a = 



atj 3 -i 

Otj 2 

a h+1 



(4.1) 



The set S a of invariant states is given by the relations: 

aj 2 -i,aj a ,ay 2 +i > (4.2) 

and 

' '1 



N 



3N 



-"72-1 



N ■ 



3N 



-a h+1 = 1. (4.3) 



We observe that S a is a 2-simplex, i. e. a triangle whose 
vertices are given by the following parameter vectors a: 



3N 
N+2 



3N 
N-2 



(4.4) 
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In order to transform to /?-space we first determine the 
matrix L by means of the formulae (IB3II - (IBf 



L = 



A', D and E, where A is given by Eq. 14. oil and: 



N-2 
3N 



(N-2){N+1) 
2JV(JV-1) 
(N+l)(N+2) 
GN(N-l) 





T 


1 N+2 




3 


V 37V 






/(N-l)(N+2) 


(N- 


-1)(JV+1) 


V 2W(7V+1) 


2(N- 


-2VN+2) 


/(JV-l)(JV-2) 


3(N- 


-1)(N+1) 


V 6N(N+1) 



The extreme points of S" are found by applying this 
matrix to the vectors given in Eq. \AAl • Since f3 is iden- 
tically equal to 1 by the normalization condition H2.1DJ1 
we can represent points in /3-space by two coordinates 
{Piifa)- One finds that S@ is a triangle in the {(3\, P2)- 
plane with the vertices: 



.4 



B 



3(N- 




/(iV- 


l)(iV- 


2) ) 


2(JV + 




2(AM 


-l)(iV4 


-2) J 


k(N- 


f 1) 

1 \ ! 




r-l)(iV- 


F2) 



(4.5) 



(4.6) 




/2(JV-2)(JV + 2) 



(N — l)(N + l) ' V (JV-1)(JV + 1) 



•7) 




Figure l: State space structure of a system composed of two 
particles with spins ji = 1 and ji = § (N = 4). The triangle 
ABC represents the set S 13 of invariant states, while the tri- 
ang le i^S" is its image under the partial time reversal. The 
polygon AA'DE represents the set Sp pt of the PPT states. 

The image i^S 1 of S@ under the partial time reversal 
is obtained by reversing the sign of the coordinate f3\. 
Consequently, Sp pt is a polygon with the four vertices A, 



A' 



>3(N-1) / (N-l)(N-2) 



2(JV + 1)' V 2(N+l)(N + 2) 



D = 0,-4 



E = 0, 



'2(N- 1){N -2) 
{N + l)(N + 2) 

I (n + i)(iy-"i)" \ 

2(N + 2)(N -2) ) 



, (4.8) 
(4.9) 
(4.10) 



Here, A' = d 2 A is the image of A under d 2 , while D and 
E are the intersections of the lines AC and AB with the 
/32-axis, respectively. The case N — 4 is illustrated in 
Fig. Similar pictures are obtained for other values of 
N. Examples are shown in Fig. Note that the origin 
of the (J3t, /?2)-plane describes the state p — j^I ® I of 
maximal entropy. 



N=4 




N=8 




N=16 




-2-1 1 -2-1 1 -2-1 1 
Pi Pi Pi 



Figure 2: The sets of the invariant states S 1 ^ and of the in- 
variant PPT states t for three different values of TV. 

To construct the set S^ ep of separable states we have 
to investigate the functions: 



(2)1 



(4.11) 



= ^E^ (1) i T i ( 9 1 V (1) )v (2) in g 2)t i/ 

8=-l 



and 



/3 2 b (1) ,^ (2) ] (4.12) 

=Vx^^ (1)|T ^ l ^ (1)) ^ (2)|T ^ )t ^ (2)) - 

g=-2 

We distinguish two cases, namely the cases of odd and of 
even N. 
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Theorem 1 For integer spins j 2 = 1,2,3,... one has 
S ppt = S@ . Hence, for all3&)N systems with odd N the 



PPT criterion represents a necessary and sufficient con- 
dition for the separability of rotationally invariant states. 

To proof this theorem we show that the vertices A, A' , D 
and E of the polygon S ppt belong to W 13 . The statement 

S@ then follows immediately from Eq. i|3.9JI . 
The point A corresponds to the parameter vector a 
given by Eq. Il3.1fi|) , It follows that this point as well as 
the point A' = # 2 A belong to W@ . Hence, it suffices to 
verify that D,E e W 3 . 

To show that E € W 13 we choose the states 



G0 



|y>W) = |l,mi = 0>, |^ 2; ) = |i2,m 2 = 0). (4.13) 

According to the selection rules for the matrix elements 
of the tensor operators l|A3ll and to Eq. l|A8ll we have 



that (<pW\T$\<p<» 



for (? = 0, ±1 and, therefore, 
01 = 0- (4.14) 



On the other hand, the non- vanishing matrix elements of 
the second-rank tensors are given by [see Eq. (IA 10ft ] : 



and 

(^ (2) |T 2 ( „ 2 V (2) > = 



2V5J2 02 + 1) 



(4.15) 



which yields: 



02 = 



y/{N + 2){N + l)N(N - 1){N - 2) 

(4.16) 



^(<P^\<P {1 Mv^\<P W ) 




(N+l)(N-l) 
2{N +2){N -2)' 



(4.17) 



We see from Eqs. SHty . 6JH and glfijl that 0\,fh) = 
E and, hence, that the point E belongs to W@. 

To show that also D belongs to W 13 we take the states 



\9 



1,0), 



\^ 2) ) = \j2,+h). 



(4.18) 



Since the state \<p^) is the same as before, Eqs. <I4 . 1 4|) 
and H4.15J1 hold true. Instead of Eq. H4.16I) . however, we 
get 



(^)|Ti 2 V (2) > = 
This gives 

h = 



2\/5[3jf - j 2 (j 2 + 1)] 



y/(N + 2)(N+ l)N(N - 1)(N - 2) 

(4.19) 



fo< 1 >|T&VW^ a V (9) > 



I2(N- 1)(N- 2) 
(A^+ l)(JV + 2) ' 



(4.20) 



A comparison with Eq. (14 . fjf) shows that (/3i,/3 2 ) = D S 
W@. This concludes the proof of the theorem. 

Let us now turn to the case of half- integer spins j 2 , 
i. e., we assume that TV is even. Of course, we again have 
that A and A' belong to W 13 . But also D e W 13 because 
the state |j 2 ,+j 2 ) exists for integer as well as for half- 
integer spins j 2 . The argument following Eq. <I4 . 1 S|> can 
thus also be applied in the present case. It follows that 
Sg ep contains at least the triangle AA'D (see Fig.|HJ). 



0.5 



-0.5 







E 


h 






I 

1 




\ 

\ y 








ppt 










D 







-0.5 





01 



0.5 



Figure 3: The set of PPT states S£ pt for N = 4. The set 
Sscp lies entirely below the straight line h through F which 
is parallel to the /3i-axis. The broken line shows the curve 
defined by Eqs. <Q7t and <Q5l . 

On the other hand, the state |j 2 ,77i 2 = 0) exists, of 
course, only for integer spins j 2 . Instead of H4.13I) we 
consider the states 



|*> (1) > = |i>o>, 

which lead to 



(2)\ 



b2,+i/2), 



0i = 0, 



02 



I {N + 2)(N-2) 
2(N+1){N- 1)' 



This shows that the point 



F= 0, 




(iV + 2)(iV-2) 
2(N+l)(N- 1) 



(4.21) 



(4.22) 



(4.23) 



belongs to W 13 . Hence, S@ contains at least the polygon 
with the vertices A, A 1 , D and F. 

We introduce the straight line h which intersects the 
point F and which is parallel to the /3i-axis (see Fig.^J. 
We are going to demonstrate that S{? ep lies entirely be- 
low this line. The line h is thus tangential to S@ ep 
and corresponds to an optimal entanglement witness (see 
Sec.0. To show this we employ the rotational invariance 
of the functions /3k [see Eq. ij3,15|) ] to obtain a suitable 
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parameterization of the states of the first spin ji = 1. 
Namely, by an appropriate rotation R any state of this 
spin can be brought into the following form: 



W) = ^|i,+i) + VT^|i,-i), 



(4.24) 



where we omit an irrelevant overall phase factor and r is a 
real parameter taken from the interval [0,1]. Invoking the 
rotational invariance we may assume without restriction 
that l^*- 1 - 1 ) is of this form. The state space of the first spin 
ji has thus only a single relevant parameter r <E [0, 1]. 

By use of the representation 14.241 the quantities Pi 
and fa become functions of the parameter r and of the 
state vector \ip^) of the second spin. Inserting Eq. 114.2411 
into Eq. I|4.11l) and using Eqs. llA8|l and ijAOp of Ap- 
pendix El we get 



fa i r , V 



(2)1 _ 



|(2r-l)(^)|7^V 2) }. (4-25) 



The function fa is found by substituting the expression 
gilH into Eq. l(4~T2l) and by using Eqs. ljAT7i)) - llA"T2"l) . 

One finds that fa can be written as the expectation value 



fa[r,^} = (^\H(X)\^) 
of the Hermitian (N x N) matrix 

H(X) = H +XH V 
Here, we have defined 



(4.26) 



(4.27) 



H ° ~ V 10 Hl 
and introduced the parameter 



T. 



(2) 



A = 2^(l-r), 



22 



< A < 1. 



T, 



(2)t 



22 



(4.28) 



For a given value of A the function fa defined by 
Eq. y.2("it is certainly smaller than or equal to the largest 
eigenvalue of H(X) which we denote by £o(A). We are 
going to demonstrate below that £o(A) is a monotoni- 
cally increasing function of A and attains its maximum 
at A = 1: 



eo(l) 



' (N + 2)(N -2) 
2(N+l)(N- 1)' 



Hence, we have 



fa[r^ {2) \ <£o(l) 



for all r and l^ 2 -*)- Note that £o(l) is equal to the /in- 
coordinate of the point F [see Eq. I|4.23|) ]. This shows 
that, as claimed, all points of W 13 and, hence, all points 
of S^ cp lie below the line h. 

To prove that £q(A) is a monotonically increasing func- 
tion of A we denote the eigenvalues of H(X) by e n (X), 




Figure 4: The largest eigenvalue £q(X) of the matrix -ff(A) 
defined by Eq. I|4.27|l for different values of TV. 



where n = 0, 1, 2 . . ., and n = labels the largest eigen- 
value. With the help of Eq. (IA6jl one verifies that H (A) 
is invariant under time reversal. It follows that if \ip) is 
an eigenstate of H(X) then also the time reversed state 
r\(p) is an eigenstate with the same eigenvalue. Since j% 
is half- integer valued the states \ip) and r\<p) are orthog- 
onal. In fact, using the antiunitarity of r and Eq. (13 . 1 If) 
we get 

(T(p\<p) = (t 2 lp\t^)* = (-if 02 {rip\ip) = -{T<p\<p), 

which shows that {np\ip) = 0. 

All eigenvalues £„(A) are thus two-fold degenerate 
and we write the corresponding eigenstates as |<^n,ifc(A)), 
where the index k = 1,2 labels the eigenstates corre- 
sponding to the same eigenvalue: |(/?„,2(A)) = r|(/? raj i(A)). 
We remark that the two-fold degeneracy is analogous to 
the Kramers degeneracy according to which the energy 
levels of an invariant system of an odd number of s pin- jr 
particles are at least two-fold degenerate (see, e. g. 

The Hellman-Feynman theorem now yields 



^ = ^ ,i(A)|^xko,i(A)). 



(4 29) I n particular, we have 



dep 
dX 



A=0 



(4.30) Q n differentiating Eq. II4.31|I once again we find: 



d 2 e = 2 IW(A)|gi|^o,i(A))| 2 > Q 



dX 



n^0,fe 



£o(A) - £„(A) 



(4.31) 



(4.32) 



(4.33) 



This shows that £o(A) is a convex function of A with 
zero derivative at A = 0. It follows that £q(A) increases 
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monotonically. Some examples of the behaviour of this 
functions are shown in Fig. 0] 

It remains to verify Eq. 114. 29J) . We first n ote t hat H (l) 
can be written with the help of Eqs. ljA10|l and ljA12j) in 

terms of the spin operator j' 2 ) 



as: 



HQ.) = 2. 



(N + 2)(N + 1)(N - 1)(N - 2) 



;(2) 



-3 



2(2) 

J y 



(4.34) 



The largest eigenvalue of this matrix is given by 



£o(l) = 2, 



(N + 2){N + l)(N - 1)(N - 2) 
3 



x 32(32 + 1) 



(4.35) 



Using N — 2j2 + 1 one shows that this equation coincides 
with Eq. g2HJ). 

We finally demonstrate that the boundary of S^ ep is 
differentiable at the point F [see Eq. Ij4.23l) ]. To this end, 
we construct a smooth curve which belongs to W^ 3 and 
passes the point F. Consider the following fixed state of 
the second spin: 



\<P 



(2)\ 



1 

V2 



if = +1/2 



v/2 



;-( 2 ) = 



-1/2). (4.36) 



This is an eigenstate of the matrix H(l) [Eq. lj4.34H l cor- 
responding to the largest eigenvalue £o(l)- Since is 
fixed the functions (3\ and depend only on the param- 
eter r and describe a curve in the (J3t, /^-plane. Writing 
r = (1 +/i)/2 and determining the matrix elements one 
finds: 



01 



3N 2 



8(N + 1)(N - 1)^ 



(l + 3/ 



(4.37) 
(4.38) 



where —1 < fx < +1. The curve described by these 
equations represents the upper half of an ellipse in the 
(Pi, P2) -plane (see Fig.EJ. It intersects the point F and 
lies entirely in W 13 . Since F is the only point of h belong- 
ing to W@, it follows that the boundary of the separa- 
bility region must be curved and that it is differentiable 
at the extreme point F, the line h being the tangent. 
Summarizing, we have shown: 



Theorem 2 For half-integer spins ji 



the 



set S£ ep of separable states is a true subset of the set of 
PPT states. Hence, for all 3 £g> N systems with even N 
the PPT criterion is only necessary and there always exist 
entangled PPT states. The line h represents the tangent 
to Ssey a t extreme point F. The set S@ ep is bounded 
by the straight lines AD and A'D and by a concave curve 
which passes the points A, A' and F. 



V. DISCUSSION AND CONCLUSIONS 

The state space structure of rotationally invariant spin 
systems has been analyzed in this paper. The set of in- 
variant states has been represented by means of two sys- 
tems of invariant operators, namely by the projections 
Pj onto the total angular momentum manifolds and by 
the invariant operators Qk composed of the spherical 
tensors. The transformation between both representa- 
tions was found to be given by a matrix L which is de- 
termined by certain 6-j symbols of Wigner. The Qk- 
representation is particularly useful in applying the PPT 
criterion for separability because the Qk are eigenoper- 
ators of the partial time reversal. The method has been 
demonstrated to lead to a complete classification of sep- 
arability of 3(g) N systems. We have shown that the PPT 
criterion is necessary and sufficient for all system with 
odd N, while entangled PPT states exist for systems with 
even N. 

Some remarks on the structure of the state space in the 
limit N — * 00 might be of interest. In this limit the value 
of the second spin j 2 becomes arbitrary large. We infer 
from Eqs. Ij4.fill - lj4.9jl that the point B then converges to 
the point A' , and C to D. At the same time F converges 
to E [see Eqs. Ij4.1fl|l and Ij4.23|l ]. Hence, as N increases 
the set Sp pt approaches the set S@ and S^ ep approaches 

S ppt . This behaviour is also indicated in Fig. [2j The 
limit N — > 00 thus corresponds to a kind of classical 
limit in which all invariant states have a positive partial 
transpose and are separable. 

The line h constructed in Sec. |^ leads to an entan- 
glement witness which we denote by W. An entan- 
glement witness is a Hermitian operator which satisfies 
tr{VVa} > for any separable state a, and trjlVp} < 
for at least one non-separable state p 0, 0|. The hy- 
perplane h corresponding to an entanglement witness W 
is defined by tr{Wp} = 0. In the case of 3 ® N sys- 
tems h is a one-dimensional line and the witness is, in 
fact, optimal 0] because h is tangential to the region 
of separable states. We have formulated the witness in 
/3-space. Transforming back to a-space one easily shows 
that the entanglement witness corresponding to h may 
be written in terms of the projections Pj as: 



W 



1 



N 



-P; 2 -l + Pj. 



1 



iV + 2' 



-P 



h + l- 



(5.1) 



This expression leads to the following physical interpre- 
tation of W. Suppose one carries out a measurement of 
the total angular momentum J on some invariant state 
p. If p is separable the inequality 



P32 



N 



P32 



P32+1 
N + 2 



> 



(5.2) 



must be satisfied, where pj = tr{Pjp} denotes the prob- 
ability of finding the value J. In other words, if the 
inequality Ij5.2jl is violated the state p must necessarily 
be entangled. 
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We exploit the witness 15.lt to design a prescription for 
the detection of entangled PPT states in 3 <g> N systems 
with even TV (bound entanglement A given state p 

is positive under partial transposition if and only if the 
corresponding point (/3i , /?2 ) lies below the line through 
A' and E, and above the line through A' and D (see 
Fig- EI- If we transform to a-space this yields the condi- 
tions 



2p* 



(iV 2 - 5)p h 



2ft 



2 + 1 



N — 1 (N + 1)(N- 1) N + l 



> 



and 



2p 



2^ 

(JV-l)(JV-2) N-l 



'h-i 



■Ph+i > °- 



(5.3) 



(5.4) 



These inequalities are equivalent to the PPT condition 
(I3.12p . Hence, entangled PPT states can be detected in 
the following way: Suppose again that a total angular 
momentum measurement is performed on some state p. 
If one finds that the measurement outcomes, i. e. the 
probabilities pj, satisfy the inequalities 115. 3J1 and 115. 4|l 
and violate the inequality i!5,2|) then the state p must be 
an entangled PPT state. 

The witness W defined in Eq. 115. 1J) does not detect all 
entangled PPT states. As has been shown in Sec. IIVI a 
part of the boundary of the region of the separable states 
is curved and, therefore, one needs an infinite number of 
linear entanglement witnesses. The upper boundary of 
S^ ep can of course be described by means of a suitable 
nonlinear equation. A possible way to derive the latter 
is to construct the envelope of appropriate families of 
curves of the type given by Eqs. H4.37J1 and H4.38I) . 

The considerations of Sec. IIVI reveal that for 3 <g> N 2 
systems half-integer spins are crucial for the emergence 
of entangled PPT states. The entanglement structure of 
systems involving half-integer spins is thus quite different 
from those with integer spins. It seems that this is closely 
connected to the fact that pure states which are invariant 
under time reversal only exist for integer spins, while for 
half-integer spins a given pure state is always orthogonal 
to its time reversed state. A clear physical interpretation 
of this result and its generalization to arbitrary N% <S> N 2 
systems is of great interest. The next step to further 
investigate this point could be to study 4 <g) N 2 systems, 
which is possible by the method developed in this paper. 
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Appendix A: SPHERICAL TENSOR OPERATORS 

We define here the irreducible spherical tensor oper- 
ators Tjc q acting on the state space of a particle 



with spin j, where N = 2j + 1, K = 0, 1, ... , 2j, and 
q = — K, . . . , +K. The tensor operators for i = 1, 2 

used in the main text are obtained by setting j = ji or 
j =h- 

The spherical tensor operators Tx q represent a com- 
plete system of operators on C N . This means that any 
operator on the state space of the spin-j particle may 
be written as a unique linear combination of the Tx q . 
Moreover, the tensors are orthonormal with respect to 
the Hilbert-Schmidt inner product: 



tr 



{4v T **} 



= 5 



KK'Oqq'- 



(Al) 



For a fixed K the (2K + 1) operators Tx q represent the 
spherical components of a tensor of rank K. They trans- 
form according to an irreducible representation of SO (3) 
which corresponds to the angular momentum K: 



+K 



dW(R)T K9 D&(R? = J2 



DfhR)T Kq - 



(A2) 



q' = -K 



For instance, the T\ q behave as components of a vector, 
and the T 2q as components of a second-rank tensor. 

The matrix elements of the tensors may be defined in 
term of Wigner's 3-j symbols as [lEHO] 



j j K 



(j,m\T K q\j,m') = V2K+l(-iy 



m —to —q 



(A3) 

The 3-j symbols are closely related to the Clebsch- 
Gordan coefficients: 



(ji,mi;j 2 ,m 2 \JM) = 
V2J + l(-l) h - j2+M 



Jl J2 
TO 2 TO2 



J 

-M 



(A4) 



According to the selection rules for the 3-j symbols the 
matrix element l|A3ll is equal to zero for to — ml — q ^ 0. 
In particular, we have Too = 

The matrix elements l|A3ll of the tensor operators are 
real and one has TjL. = Tj£ g = {—l) q TK,- q . It fol- 
lows that the Tn q are eigenoperators of the time reversal 
transformation § which was defined in Eq. IjS.lOp . In 
fact, using the transformation behaviour i|A2j) of the ten- 
sors and the fact that a rotation by 7r about the y-axis is 
represented by the unitary matrix 



4?w 



-1) 



K-q' 



finds 



dT Kq = VTl W = {-l) K T t 



Kg- 



(A5) 



(A6) 



As a consequence the operators Qk which have been in- 
troduced in Eq. 112. 7J1 are eigenoperators of the partial 
time reversal $2 = I ® 



$2Qk = (-1) k Qk 



(A7) 
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We finally list the non-vanishing matrix elements of 
the tensor operators needed in Sec. IIVI 



(j,m\T 10 \j,m) = 2m, 



(j,m\ T n\j> m + 1 ) = -\ 



N(N — 1)(JV + 1)' 



' 6(j — m)(j + m + 1) 
N(N -1)(N +1) 



(A8) 



(A9) 



(j,m\T 20 \j,m) 



2V5[3m 2 -j(i + l)] 



\J (N + 2)(N + 1)N(N - 1)(JV - 2) ' 

(A10) 



(j,H T 2lli) TO + 1) = 
-^5(1 + 2m) 



(All) 



6(j — m)(j +771 + 1) 



(JV + 2) (AT + 1)N(N - 1)(JV - 2) ' 



<J, "T-I^^lj, 777. + 2) 



(A12) 



V5< 



'6(j — 771 — l)(j — m)(j + m + l)(j + m + 2) 
(AT + 2) (AT + 1)AT(AT - l)(N - 2) ' 



Appendix B: PROOF OF RELATION l(234|) 

The starting point is given by Eq. 12.13L We insert 
into this equation the definitions <2.2ll and l|2.7|i for the 
invariant operators Pj and Qk , and introduce complete 
sets of product basis states |7i,mi; 7 2 ,to 2 ). This yields a 
multiple sum over products of two Clebsch-Gordan coef- 
ficients and two matrix elements of the tensor operators. 
By use of Eqs. l!A3j) and QA4I) the Clebsch-Gordan co- 
efficients as well as the matrix elements of the spherical 
tensors can be written in terms of the 3-j symbols. We 



also use the selection rules for the 3-j symbols and their 
symmetry properties. This procedure leads to the follow- 
ing sum over 4-fold products of 3-j symbols: 



Lkj = 



y/(2K + 1)(2 J + l)(-l) h+j3+J x 
X({mi}) x 



(Bl) 



jl 32 J 
mi 777,2 ^3 

32 K 

-7722 777 6 



32 
-777-4 



Jl Jl 
-777l ?775 

.72 
7774 



K 

m 6 
31 J 

-777 5 -777 3 



where x({ m i}) ls a phase factor: 

x({"*i}) = (-l)^+ mi (-lp +m2 (-l) J+m3 x 

(_\\h+mi f_-j\ji+me f_-\\K+me 

The sum over the quantum numbers mi, ...,iri6 in 
Eq. JB II) exactly corresponds to a certain 6-j symbol of 
Wigner [^j. A general 6-j symbol involves six angular 
momenta and is written as 



Jl J2 J3 

h h J6 



(B2) 



The sum of Eq. I|B1J) is equal to the 6-j symbol l|B2j) with 
js = J, ji = j2, 35 = ji and je = K. Hence, we see that 
Eq. IB 11 reduces to Eq. 12.141 . We remark that a similar 
technique has been used in Ref. jjjjj m order to derive 
an expression for the matrix which represents the partial 
time reversal $2 in the Pj-representation. 

By use of the formulae for the 6-j symbols j2^| we find 
that the first three rows of L are given by 



2J+ 1 

AiA 2 : 



(B3) 



and 



Lu = -2^/3(27+1) 



Jl(jl+1)+J2(j2 + 1)- J(J+l) 

y/(N! - l)iVi(JVi + 1)(JV 2 - 1)N 2 (N 2 + I)' 



(B4) 



L 2J = 2^5(2J+1) 



3X(A-l)-4jl(7l + l) 7 2 (72 + l) 



y/(N! - 2)(JVi - l)Ai(A^i + l)(JVi + 2){N 2 - 2){N 2 - l)iV 2 (JV 2 + 1)(JV 2 + 2) ' 



(B5) 



where X = ji(ji + 1) + j 2 (j 2 + 1) - J {J + 1). 
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